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Q_^. Abstract 

2 i The effect of massive neutrinos on matter power spectrum is discussed in the context of f{R) 

gravity. It is shown that the anomalous growth of density fluctuations on small scales due to the 

scalaron force can be compensated by free streaming of neutrinos. As a result, models which predict 
(N ; 
^ I observable deviation of the equation-of-state parameter wbe from wde = — 1 can be reconciled with 
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observations of matter clustering if the total neutrino mass is 0(0.5 eV). 



f{R) gravity is a simple and nontrivial extension of General Relativity (GR) which has 
recently received much attention as an alternative to the A-Cold-Dark-Matter (ACDM) 
model that can explain the present cosmic acceleration. The basic idea is to modify the 
Einstein-Hilbert action by using the general function of Ricci scalar R as 

s = Jd''x^rgf{R) + s^, (1) 

where Sm indicates a matter action (see Ref. [1] for a recent review and an extensive list 
of references). If f"{R) is not zero identically, this modified gravity contains an additional 
degree of freedom corresponding to a new scalar field, dubbed scalaron in Ref. [2] where 
the simplest, f{R) = R + -R^/6M^, variant of such theory (with small one-loop quantum 
gravitational corrections) was used to construct an internally self-consistent inflationary 
model of the early Universe with a graceful exit from inflation to the subsequent radiation- 
dominated Friedmann-Robertson- Walker (FRW) stage after the gravitational creation of 
matter and reheating. However, f{R) gravity can be used to describe dark energy in the 
present Universe as well, although with the use of more complicated functions f{R) having 
a nontrivial structure for low values of R. 

To have the correct Newtonian limit for R ^ Rq = R{to) ~ H^, where to is the present 
moment and Hq is the Hubble constant, as well as the standard matter-dominated FRW 
stage with the scale factor behaviour a{t) oc t"^^^ driven by cold dark matter and baryons, 
the following conditions should be fulfilled: 

i/(i?) - i?i « i?, i/'(i?)-i|«i, i?r(i?)«i, i?»i?o, (2) 

where the prime denotes a derivative with respect to the argument R. The same conditions 
guarantee smallness of non-GR corrections to a space-time metric for more general back- 
grounds for which GR has to be used in full. In addition, the stability condition f"{R) > 
has to be satisfied to guarantee that the standard matter- and radiation-dominated FRW 
stages remain attractors with respect to an open set of neighboring isotropic cosmological 
solutions in f{R) gravity. In quantum language, this condition means that the scalaron is 
not a tachyon. Note that the other stability condition, f'{R) > 0, which means that gravity 
is attractive and the graviton is not a ghost, is automatically fulfilled in this regime. Spe- 
cific functional forms that satisfy these and other necessary conditions have been proposed 
in Refs. [3-5], and much work has been conducted on their cosmological consequences. 



K - \kr = -svrG ( t;;,. + t,%^. ) , (3) 



We can express field equations derived from action (1) in tlie following Einsteinian form, 

where 

SttGT^^j.e) ^{F- l)i?^ - i(/ - R)6i: + (V^V, - 5^,n)F, (4) 

and we have also defined F{R) = f'{R) = df /dR. Working in the spatially fiat Friedmann- 
Robertson- Walker (FRW) space-time with the scale factor a{t), we find 

3H^ = 8nGp - 3(F - 1)H^ + ]-{FR - f) - 3HF, (5) 

2H = -SnGp - 2{F - l)H - F + HF, (6) 

where H is the Hubble parameter and p is the energy density of the material content, which 
we assume to consist of nonrelativistic matter. From these expressions, the effective DE 
equation-of-state parameter wbe = -Pde/pde is given by 

-2H - MGp 



1 + wde 

(7) 



3H^ - SttGp 
2(1 - F){-d/a + H^) + F-HF 



3(1 - F){--d/a) + {R- /)/2 - 3HF 
Note that the phantom crossing is naturally realized in viable f{R) theories as shown in 
Refs. [3, 6-8]. 

We specifically adopt the following functional form that satisfies all the above 
requirements [5], 



f{R) =R + \Rs 



-|V^. 



(8) 



where n. A, and Rs are model parameters. It is sufficient to describe the recent evolution 
of the Universe, although it has to be modified both for R < Rq (including the region 
R < 0) and for very large positive R to avoid problems in the early Universe cosmology (see 
Refs. [6, 9]). 

In f{R) gravity, the evolution equation for density fluctuation with comoving wave num- 
ber k in Fourier space, (5k, in the deeply subhorizon regime is given by [10, 11] 

5k + 2i7(5k - 47rGeffp5k = 0, (9) 

where 



Gl + 4_|^ 

'a2 F 



G.. - ^ytM- ^'°^ 



In the high-curvature regime, the differential equation (9) can be solved analytically [9], 
and we analyzed it in the general regime in terms of numerical calculation in the previous 
paper [8]. It is revealed that the effective gravitational coefficient Geff enhances the growth 
of density fluctuations compared with that in the standard ACDM model due to an extra 
force mediated by the scalaron. To suppress the deviation of the matter power spectrum 
from that in the ACDM model within, say 10%, the parameter space for n and A should be 
restricted. For example, A should be greater than 8 in n = 2 [8]. 

In this paper, we show that the anomalous growth of density fluctuations on small scales 
due to the scalaron is rectified if neutrino masses are taken into consideration. That is, the 
free streaming of neutrinos partially erases small-scale density fluctuations, which can be 
harmful in Einstein gravity, to compensate their unwanted extra growth in f{R) gravity. 

To set up the initial condition, we calculate the evolution of density fluctuations for 
z > 10 using CAMB [12, 13], and obtained P{k,z = 10), power spectrum aX z = 10. We 
assume one mass eigenstate for neutrino, and all three neutrinos have practically the same 
masses. Note that the density parameter of neutrinos is 

f 1.685 X 10-^ (^m^ = 0) 
n,ec^{ y^ (11) 

'^ 94.1 eV ^^ ^^ ^' 

and the free streaming of massive neutrinos suppresses fluctuations below the scale [14, 15] 

U^) - T^^ (^ (^] ''' h Mpc-. (12) 

^ ' (l + z)i/2 Vl eVJ \0.27J ^ ^ ' 

Therefore, (5k in Eq. (9) implies density fluctuations for CDM below that scale. 

At 2; > 10, the effect of modifled gravity is sufficiently small. Therefore, we start to solve 
the evolution of density fluctuations in f{R) gravity using Eq. (9) at 2; = 10. We have found 
that the physical wave number crossing the scalaron mass today is given by 

^ ^ 3.2 Ai-««/fo, (13) 

for n = 2[8], which is to be compared with Eq. (12). It takes 1.07 x lO^'^/i and 8.44 x 
10~^h Mpc~^ for A = 1 and 3, respectively. 

The present power spectrum is constructed as 

P{k, z = 0) = P{k,z = 10) i^j^^^'j . (14) 
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FIG. 1: Deviation of power spectrum for n = 2, A = 1, 3. Captions indicate total neutrino mass, 
and it is assumed that three neutrinos have the same masses. 
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FIG. 2: Left: Relative deviation of linear matter power spectrum from ACDM model at fc = 
0.174/i Mpc^ as a function of A for various values of neutrino mass. Right: The region between 
two lines for each n corresponds to the parameter space where the deviation of power spectrum is 
smaller than 10%. 

Figure 1 shows cancellation of the contributions from j{K) gravity and neutrino masses. 
Equations (12) and (13) explain the threshold wave number. In particular, the transition 
occurs at a similar wave number for A = 3 when we take n = 2. The deviation at A; = 
0.174/i Mpc~^, which is the scale corresponding to erg normalization [8], is minimum when 
the total neutrino masses are 0.6 and 0.5 eV for A = 1 and 3, respectively. 

In the left panel of Fig. 2, the relative deviation of the power spectrum from that of the 
ACDM model is depicted as a function of A for various total neutrino masses with n = 2. The 
right panel of Fig. 2 shows the range of total neutrino masses where the relative deviation 
of the power spectrum at k = 0.174/i Mpc~^ remains smaller than 10%. 

As seen in the right panel, models with smaller values of A, which tend to exhibit a 
larger deviation from the ACDM model, become compatible with the observation of flue- 
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FIG. 3: Possible range of time variation of u;de (left) and growth index 7 (right) as a function of 
total neutrino mass. 7 is measured at the comoving wave number k = 0.174/i Mpc~ . 

tuations. In other words, models with deviations of wbe and the growth index 7(2;) = 
log ( 6/H6 ] I log fim large enough to be observable become viable thanks to the suppression 
of small-scale fluctuations due to finite neutrino masses. Figure 3 depicts the maximum 
range of the time variation of wde and 7 that can be realized for each value of the total 
neutrino mass. That is, if we choose a minimum possible value of A for each total neutrino 
mass, wde varies between the upper and lower edges of the spatula shape at the correspond- 
ing neutrino mass. Similarly, 7 evolves between the upper and lower edges of the ax-shape. 
Note that 7 takes an almost constant value, 0.55, in the ACDM model, which corresponds 
to the upper bound in f{K) gravity. 

Finally, we comment on how nonlinear effects may modify our results. So far, nonlinear 
effects on matter clustering in /(-R) gravity have been analyzed in higher order perturbation 
theory [16] and N-body simulations [17]. The latter has shown that nonlinear effects lower 
the relative deviations of the power spectrum, (PfRc — -Pacdm)/-Pacdm, by about 5% in a 
specific case f{K) = i? — 2A — fR^B^/R with |//jo| = 10^'*, where _Ro is the current value of 
the scalar curvature. It has also been shown that the nonlinear effects tend to suppress the 
anomalous growth of small-scale fluctuations observed in linear perturbation in f{R) theory. 
We may therefore conclude that nonlinear effects will not change our results qualitatively. 

In conclusion, neutrino rest-masses with ^^Tn^, = 0(0.5 eV) relax the most critical 
constraint on f{R) gravity, which follows from the anomalous /c-dependent growth of density 
perturbations in the cold dark matter + baryon component at recent redshifts, making 
possible for wde and '~f{k) deviate from those of the ACDM model noticeably. One can 



distinguish the effects of f{R) and massive neutrinos by analyzing the ratio of correlation 
between galaxies and curvature perturbation probed by weak lensing to that between galaxies 
and velocity field, namely, an estimator Eq = VLm/{FP)^ where (5 is the growth rate[18]. 
According to the current observational analysis, both Einstein gravity and f{R) gravity are 
consistent with the data. However, in the future, one can break the degeneracy in principle. 
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